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We discuss recently generated high-temperature series expansions for the free 
energy and the susceptibility of random-bond q-state Potts models on hypercu- 
bic lattices. Using the star-graph expansion technique quenched disorder aver- 
ages can be calculated exactly for arbitrary uncorrelated coupling distributions 
while keeping the disorder strength p as well as the dimension d as symbolic 
parameters. We present analyses of the new series for the susceptibility of the 
Ising (q = 2) and 4-state Potts model in three dimensions up to order 19 and 
18, respectively, and compare our findings with results from field-theoretical 
renormalization group studies and Monte Carlo simulations. 
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1. Introduction 

Systematic series expansions [1] for statistical physics models defined on a lat- 
tice provide an useful complement to field-theoretical renormalization group studies 
and large-scale numerical Monte Carlo simulations. This is in particular true when 
studying phase transitions and critical phenomena of quenched, disordered systems. 
In the field-theoretic treatment [2] the necessary average over disorder realizations 
at the level of the free energy requires the application of the so-called "replica trick" 
which loosely speaking introduces n different, interacting copies of the original sys- 
tem, with the formal limit n — > taken at the end. In the numerical approach the 
average over a large but finite number of different disorder realizations can, at least 
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in principle, be performed explicitly but is very time consuming such that only few 
points in the vast parameter space of the systems can be sampled with realistic 
effort. Moreover extrapolations of the data on finite lattices to the infinite- volume 
limit are required. Using high-temperature series expansions, on the other hand, 
one can obtain for many quantities exact results up to a certain order in the inverse 
temperature. Here the quenched disorder is treated exactly and the infinite-volume 
limit is implicitly implied. Moreover, one can keep the disorder strength p as well 
as the dimension d as symbolic parameters and therefore analyse large regions of 
the parameter space of disordered systems. The critical part of the series expan- 
sion approach lies in the extrapolation techniques which are used in order to obtain 
information on the phase transition behaviour from the finite number of known co- 
efficients of the high-temperature series. While for pure systems this usually works 
quite well, one can question the use of these extrapolation techniques in disordered 
systems, where the singularity structure of the free energy or susceptibility may be 
very complicated, involving Griffiths-type singularities or logarithmic corrections [3]. 

Pure Potts models show either first- or second-order phase transitions, depending 
on the dimension d and the number of states q. Since in the second-order case the 
specific-heat exponent a is non-negative for this class of models, the Harris criterion 
[4] suggests for the corresponding disordered systems either the appearance of a new 
random fixed point (d = 2, q = 3,4 and d = 3, q = 2) or logarithmic corrections 
to the pure fixed point (d — 2, q — 2). At first-order phase transitions, randomness 
softens the transitions [5] . For d = 2 even infinitesimal disorder induces a continuous 
transition [6,7], whereas for d = 3, q > 2 a tricritical point at a finite disorder 
strength is expected [8]. 

In this work we studied these scenarios by means of "star-graph" high-temper- 
ature series expansions where the disorder average can be taken at the level of indi- 
vidual graphs. Using optimized cluster algorithms for the symbolic, exact calculation 
of spin-spin correlators on finite graphs with arbitrary inhomogeneous couplings, we 
obtained series expansions for the free energy and susceptibility in the inverse tem- 
perature up to order 19 respectively 18 for bond-diluted Ising and Potts models in 
dimensions d < 5, and up to order 17 in arbitrary dimensions. Here we shall fo- 
cus on analyses of these series in three dimensions where a direct comparison with 
field-theoretic renormalization group studies and recent Monte Carlo simulations is 
possible. 



The ferromagnetic disordered g-state Potts model on hypercubic lattices Z rf is 
defined by the partition function 
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where (3 = 1/ksT is the inverse temperature, Jij are quenched (non-negative) 
nearest-neighbour coupling constants, the spins can take the values Sj = 1, . . .,q, 
and 8.. is the Kronecker symbol. In our series expansion the combination 



e f3Jij _ i 
e pJij — i _|_ g 



(2) 



will be the relevant expansion parameter which in the Ising case (q = 2) simplifies 
to = tanh(/3Jjj/2). In the symmetric high-temperature phase, the susceptibility 
associated with the coupling £\ hi(q8 Sit i — — 1) to an external field hi is given 
for a graph with N spins by summing over all two-point correlations, 
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Here the brackets 



indicate the quenched disorder average which in our case 



is taken over an uncorrected bimodal distribution of the form 



P(J l3 ) = (l-p)S(J ij - Jo) +pS(Jij ~ RJo). (4) 

Besides bond dilution (R = 0), which will be in the focus of the present work, 
this also includes random-bond ferromagnets (0 < R < 1) and the physically very 
different class of spin glasses (R — —1) as special cases. Other distributions such as 
Gaussian distributions can, in principle, also be considered with our method. 



3. Series generation methodology 

In this section we briefly review the main technical ingredients necessary for 
our high-temperature series study. We begin with a few basic notations from graph 
theory. A graph of order E consists of E links connecting N vertices. We consider 
only connected, undirected graphs that are simple: no link starts and ends at the 
same vertex (no tadpoles) and two vertices are never connected by more than one 
link. Subgraphs are defined by the deletion of links. In this process, isolated vertices 
can be dropped. Since each link may be present or absent, a graph of order E 
has 2 E (not necessarily non-isomorphic) subgraphs. These subgraphs may consist of 
several connected components and are called clusters. If the deletion of one vertex 
renders the graph disconnected, such a vertex is termed articulation point. The 
"star graphs" we are considering here are thus just defined by the absence of such 
articulation points. A graph is bipartite if the vertices can be separated into red and 
black vertices so that no link connects two vertices of the same color. Equivalently, 
all closed paths in the graph consist of an even number of links. Hypercubic lattices 
are evident examples for bipartite graphs. 

There are a couple of well-established methods [1] known for the systematic 
generation of high-temperature series expansions which differ in the way relevant 
subgraphs are selected or grouped together. A recently developed alternative method 
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[9] exploits ideas from so-called finite-lattice methods, usually employed before for 
the generation of /ou»-temperature series. Using a clever reformulation of the method, 
Arisue et al. [10] succeeded to generate a very impressive 32th order world-record 
high-temperature susceptibility series for the pure Ising model in three dimensions. 

For the class of classical 0(N) spin models without disorder, quite long series 
(up to order j3 25 ) have also been produced by linked-cluster expansions [11]. This 
technique also allows one to obtain series for more involved observables (such as the 
second moment of the spin-spin correlation function yielding the correlation length) 
which have no star-graph expansion. Furthermore, it works with free embeddings 
of graphs into the lattice which can be counted orders of magnitude faster than 
the weak embedding numbers needed by the star-graph technique. Nonetheless, the 
linked-cluster method has not yet been applied to problems with quenched disorder. 

The star-graph method can be adopted to systems involving quenched disorder 
[12,13] (as also can the no-free-end method [14]) since it allows one to take the 
disorder average on the level of individual graphs. The basic idea is to assemble 
the value of some extensive thermodynamic quantity F on a large or even infinite 
graph from its values on subgraphs: Graphs constitute a partially ordered set under 
the "subgraph" relation. Therefore, for every function F{G) defined on the set of 
graphs exists another function Wf(G) such that F(G) = J2 g cc Wf{q), for all graphs 
G. This function can be calculated recursively via W F (G) = F(G) — J2 gcG W F (g), 
resulting for an infinite (e.g. hypercubic) lattice in F(Z d ) = J2 G (G : Z d )WF(G), 
where (G : Z d ) denotes the weak embedding number of the graph G in the given 
lattice structure [15]. 

The following observation makes this a useful method: Let G be a graph with an 
articulation vertex where two star subgraphs G 1>2 are glued together. Then W F (G) 
vanishes if F{G) = F(G 1 ) + F(G 2 )- An observable F for which this property is 
true on arbitrary graphs with articulation points allows a star-graph expansion. All 
non-star graphs have zero weight Wp in the sum for F(Z d ). It is easy to see that the 
(properly normalized) free energy log Z has this property and it can be proved [13] 
that the inverse susceptibility 1/x has it, too, even for arbitrary inhomogeneous 
couplings Jij. This restricts the summation for F(Z d ) to a sum over star graphs. 
The linearity of the recursion relations then enables the calculation of quenched 
averages over the coupling distribution on the level of individual graphs. 

The resulting recipe for the susceptibility series is: 

• Graph generation and embedding number counting. 

• Calculation of Z(G) and the correlation matrix 
M nm (G) = Tr (q5(S n , S m ) - l) e -^«^» 

for all graphs as polynomials in E variables % defined in (2). 

• Inversion of the Z polynomial as a series up to the desired order. 

• Averaging over quenched disorder, 
N nm {G) = [M nm /Z] p{J) , 

resulting in a matrix of polynomials in (p,v). 
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Table 1. Number of star graphs with E > 8 links and non-vanishing embedding 
numbers on Z rf . For E = 1,4, 6, and 7 only a single star graph exists. 



order E 


8 9 10 


11 


12 


13 


14 


15 


16 


17 


18 


19 


# 


2 3 8 


9 


29 


51 


142 


330 


951 


2561 


7688 


23078 



• Inversion of the matrix N nm and subgraph subtraction, 

• Collecting the results from all graphs, 
1/X = E G (G:Z-)W X (G). 

Algorithmically the most cumbersome part of this recipe is the first step, the 
generation of star graphs and calculation of their (weak) embedding numbers. The 
graph generation is usually done by recursively adding nodes and edges to a list of 
smaller graphs. To make sure that no double counting occurs this requires an isomor- 
phism test, i.e., the decision whether two given adjacency lists or adjacency matrices 
describe the same graph modulo relabelling and reordering of edges and nodes. We 
employed the NAUTY package by McKay [16] which allows very fast isomorphism tests 
by calculating a canonical representation of the automorphism group of the graphs. 
By this means, we classified for the first time all star graphs up to order 19 that 
can be embedded in hypercubic lattices, see Table 1. As with any series expansion, 
the effort grows exponentially with the maximal order of the expansion, rendering 
each new order roughly as "expensive" as all previous orders taken together. This is 
illustrated in Fig. 1 where already the number of star graphs is seen to grow expo- 
nentially as a function of the links E. The exponential fit in the range E — 13 - 19 
suggests that the number of star graphs increases roughly by a factor of 2.8 in each 
of the next higher orders, predicting about 65 000 different star graphs with E = 20 
and about 180 000 with E = 21. 

For each of these graphs we calculated their (weak) embedding numbers for d- 
dimensional hypercubic lattices (up to order 17 for arbitrary d, order 18 (general 
g-state Potts) and 19 (Ising) for dimensions d < 5). Two typical results are de- 
picted in Fig. 2. For the embedding count we implemented a refined version of the 
backtracing algorithm by Martin [15], making use of a couple of simplifications for 
bipartite hypercubic lattices Z d . After extensive tests to find the optimal algorithm 
for the "innermost" loop, the test for collisions in the embedding, we ended up using 
optimized hash tables. 

The second step of the series generation requires the exact calculation of the 
partition function and the matrix of correlations M nm for each star graph with ar- 
bitrary symbolic couplings Jy defined on the E < 19 edges. The crucial observation 



5 



Hellmund and Janke 




Figure 1. Growth behaviour of the number of star graphs with E links that can 
be embedded in hypercubic lattices 7L d . 



is that this can be done most efficiently by using the cluster representation 
Z oc Z = q~ N Tr JJ [l - % + Vijq5 SuSj ] 




" %)J , (5) 

where the sum goes over all clusters C C G, e is the number of links of the cluster and 
c the number of connected components of C. The reduced partition function Z = 
Zq E ~ N / Y\(ij)i. el3Jx3 ~ 1 + 9) i s normalized such that log Z has a star-graph expansion. 
Similarly, the calculation of the susceptibility involves the matrix of correlations 

M nm oc Q e+ °~ N [ II % I [ II (! - I » ( 6 ) 



r?. t?7. 



where the sum is restricted to all clusters C nm C G in which the vertices n and m 
are connected. 

This representation essentially reduces the summation over q N states to a sum 
over 2 E clusters which, compared with previous implementations, results is a huge 
saving factor in computing time (of the order of 10 6 ). Further improvements result 
if the 2 E clusters belonging to a graph are enumerated by Gray codes [17] such that 
two consecutive clusters in the sum (5) differ by exactly one (added or deleted) link. 
In the Ising case q = 2 another huge simplification takes place since only clusters 
where all vertices are of even degree contribute to the cluster sum. 
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7620(2) + 76851600(3) + 14650620864 (J) 
+ 404500471680(g) + 3355519311360(^) 




12048(3) + 396672 (J) + 2127360(5) + 24 88320(g) 



Figure 2. Two star graphs of order 17 and 19 and their weak embedding numbers 
up to 6 dimensions. 

Since general purpose software for symbolic manipulations turned out to be too 
slow for our purposes, we developed a C++ template library using an expanded 
degree-sparse representation of polynomials and series in many variables. For arbi- 
trary-precision arithmetics the open source library GMP was used. Finally, for the 
case of bond dilution (R = in (4)) considered here, we made use of the fact that 
the disorder average is most easily calculated via 



4. Series analysis: techniques and results 

4.1. Bond-diluted 3D Ising model 

Disordered magnetic systems belonging to the 3D Ising model universality class 
have been studied extensively in experiments [18-20] and also by field theoretical and 
numerical methods. A comprehensive compilation of recent results can be found in 
[21], showing a wide scatter in the critical exponents of different groups, presumably 
due to large crossover effects. 

Our high-temperature series expansion for the susceptibility up to order 19 is 
given with coefficients as polynomials in p, x( v ) = En^W" 11 [22]- Therefore it 
should be well-suited for the method of partial differential approximants [23] which 
was successfully used to analyse series with an anisotropy parameter describing the 
crossover between 3D Ising, XY and Heisenberg behaviour [24]. But this method 
was unable to give conclusive results. Therefore we confined ourselves to a single- 
parameter series for selected values of p. 

The ratio method assumes that the expected singularity of the form 




(7) 



X (y) = A(V C -V) 7 + . . . 



(8) 
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Figure 3. Ratio approximants for different dilutions p vs. 1/n. In order to make 
them visually comparable, they are (except for p = 0.75) normalized by their 
respective critical couplings v c . 



is the closest to the origin. Then the consecutive ratios of series coefficients behave 
asymptotically as 

r. = -!± = (l + 2zl). (9) 



-l V n 

Figure 3 shows these ratios for different values of p. For small p they show the 
typical oscillations related to the existence of an antiferromagnetic singularity at 
— v c . Near the percolation threshold at p c = 0.751 188 [25] (where T c goes to 0, v c 
to 1) the series is clearly ill-behaved, related to the exp(l/T) singularity expected 
there. Besides that, the slope (related to 7) is increasing with p. 

The widely used DLog-Pade method consists in calculating Pade approximants 
to the logarithmic derivative of x( v )i 

^M = ^- + .... (10) 
dv v c — v 

The smallest real pole of the approximant is an estimation of v c and its residue 
gives 7. The results presented in Table 2 are the averages of 45 - 55 different Pade 
approximants for each value of p, with the error in parentheses indicating the stan- 
dard deviation. The scattering of the Pade approximants increases with p, getting 
again inconclusive near the percolation threshold. Nevertheless, up to about p = 0.6 
the series estimates for v c respectively T c are in perfect agreement 1 with the Monte 
Carlo (MC) results of Ref . [26] . This is demonstrated in Fig. 4 where also the (prop- 
erly normalized) mean-field and effective-medium approximation [27] are shown for 
comparison. 

The critical exponent 7, as provided by this method, apparently varies with 
the disorder strength. More sophisticated analysis methods, such as inhomogeneous 

1 Notice that "p" in the present notation corresponds to "1 — p" in Ref. [26]. 
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Table 2. Transition points v c = tanh(/3 c Jo/2) and critical exponents 7 for different 
dilutions p as obtained from DLog-Pade approximants. 



p 


V c 


7 





0.21813(1) 


1.2493(7) 


0.075 


0.23633(1) 


1.2589(8) 


0.15 


0.25788(1) 


1.2714(8) 


0.225 


0.28382(1) 


1.2873(10) 


0.3 


0.31566(2) 


1.305(4) 


0.375 


0.35557(5) 


1.329(4) 


0.45 


0.40743(10) 


1.365(6) 


0.525 


0.4772(2) 


1.400(10) 


0.6 


0.576(1) 


1.435(60) 



differential approximants [28,29], the Baker- Hunter method [30] or the methods Ml 
and M2 [31], especially tailored to deal with confluent singularities as one would 
expect in a crossover situation, give improved results in the pure (p = 0) case but 
do not essentially change the results in the presence of disorder. 

Thus, while for theoretical reasons we still find it likely that the variation of 7 
with the disorder strength can be attributed to neglected or insufficiently treated 
correction terms, it proved clearly impossible to verify this effect in the series anal- 
ysis. In fact, a plot of 7 vs. p does not even show an indication of a plateau. In the 
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Figure 4. Transition temperatures of the bond-diluted Ising model for different di- 
lutions p as obtained from our DLog-Pade high-temperature series (HTS) analyses 
and from Monte Carlo (MC) simulations [26]. For comparison also the (properly 
normalized) mean-field and effective-medium approximations are shown. 
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central disorder regime, p = 0.3 - 0.5, the high-temperature series estimates given in 
Table 2 are at least compatible with Monte Carlo results for site and bond dilution 
[26,32,33] which cluster quite sharply around 7mc = 1-34(1). Field-theoretic renor- 
malization group estimates [21,34] favor slightly smaller exponents of 7rg = 1-32 - 
1.33, while experiments [18-20] report values between 7 cxp = 1.31 - 1.44, cp., e.g., 
the table in Ref. [35]. 

4.2. Bond-diluted 4-state Potts model 

In three dimensions the 4-state Potts model exhibits in the pure case a strong 
first-order transition [36] which is expected to stay first order up to some finite 
disorder strength, before it gets softened to a second-order transition governed by a 
disorder fixed point. 

In the latter regime we are interested in locating power-law divergences of the 
form (8) from our susceptibility series up to order 18 [37,38] . To localize a first-order 
transition point, however, a high-temperature series alone is not sufficient since there 
the correlation length remains finite and no critical singularity occurs. In analysing 
series by ratio, Pade or differential approximants, the approximant will provide an 
analytic continuation of the thermodynamic quantities beyond the transition point 
into a metastable region on a pseudo-spinodal line with a singularity T* < T c and 
effective "critical exponents" at T*. Again we first employed the ratio method which 
is the least sophisticated method of series analysis, but usually it is quite robust 
and gives a good first estimate of the series behaviour. Figure 5 shows these ratios 
for different values of p. They behave qualitatively similar to the Ising model case 
(oscillations caused by the antiferromagnetic singularity at — v c , strong influence of 
the percolation point at p c ~ 0.75). Notice that the slope (oc 7 — 1) is increasing 
with p, changing from 7 < 1 to 7 > 1 around p = 0.5. 

Figure 6 compares the critical temperature, estimated from an average of 25 - 
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1.04 
J= 1 .02 
1 

0.98 
0.96 
0.94 

0.05 0.1 0.15 0.2 

1/n 

Figure 5. Ratio approximants for different dilutions p vs. 1/n (normalized by v c 
as in Fig. 3). 




10 



High-temperature series expansions for random Potts models 



1.6 












MCdata — i — - 
series data — •*— 


1.4 












- 


1.2 
















0.03 
" 0.02 












1 












0.8 


0.01 













0.6 


" -0.01 















0.2 0.4 


0.6 





0.1 0.2 0.3 0.4 0.5 0.6 

P 



Figure 6. Transition temperatures of the bond-diluted 4-state Potts model for 
different dilution p as obtained from Monte Carlo (MC) simulations [39] and 
DLog-Pade series analyses. The inset shows the difference between the two esti- 
mates. 



30 Pade approximants for each value of p, 2 with the results of recent Monte Carlo 
simulations [39]. For small p, in the first-order region, the series underestimates the 
critical temperature. As explained above, this is an estimate not of T c but of T*. 
Between p = 0.3 and p = 0.5, the estimates confirm, within errors, the Monte Carlo 
results, indicating that now both methods see the same second-order transition. 
Beyond p — 0.5, the scatter of different Pade approximants increases rapidly, related 
to the crossover to the percolation point. 

The situation is more complicated with respect to the critical exponent 7. The 
DLog-Pade analysis gives inconclusive results due to a large scattering between 
different Pade approximants, as shown in Fig. 7. One possible reason for this failure 
is the existence of confluent singularities. The dots in Eq. (8) indicate correction 
terms which can be parametrized as follows: 

X (v) = A(v c - u)-t[1 + A 1 (v c - v) A ^ + A 2 (v c - v) A * + ...], (11) 

where A, are the confluent correction exponents. Among the various sophisticated 
analysis methods (inhomogeneous differential approximants [28,29] and the methods 
Ml and M2 [31]), in the case at hand, the Baker-Hunter method [30] appeared to be 
the most successful, giving consistent results at larger dilutions p > 0.35 where the 
leading-term DLog-Pade analysis failed. The Baker-Hunter method assumes that 
the function under investigation has confluent singularities 

N / A" Al 

n*) = j> 1-- =E fl ^> ( 12 ) 

i=l V ZcJ n=0 

2 Again, "p" in the present notation corresponds to "1 — p" in Ref. [39]. 
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Figure 7. Scattering of different Pade approximants at dilution p = 0.4: critical 
exponent 7 against critical coupling v c . 
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*- 0.965 - 
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95 / 

0.32 0.3205 0.321 0.3215 0.322 0.3225 0.323 0.3235 0.324 0.32 0.3205 0.321 0.3215 0.322 0.3225 0.323 0.3235 0.324 

Figure 8. Values for the critical exponent 7 and amplitude A at p = 0.4 as 
function of trial v c estimates from the Baker-Hunter analysis. From the clustering 
of different Pade approximants in both pictures we estimate v c = 0.3217, 7 = 
0.966, and A = 1.21. 

which can be transformed into an auxiliary function g{t) that is meromorphic and 
therefore suitable for Pade approximation. After the substitution z = z c (l — e~ l ) we 
expand F(z(t)) = J2 n c n t n and construct the new series 

N . 

g(t) = J2n\c n t n = J2 T ^ rt , (13) 

n=0 i=l 1 

such that Pade approximants to g(t) exhibit poles at t — 1/Aj with residues — Ai/\. 
This method is applied by plotting these poles and residues for different Pade ap- 
proximants to g{t) as functions of z c . The optimal set of values for the parameters 
is determined visually from the best clustering of different Pade approximants, as 
demonstrated in Fig. 8. 

Using this method, our results for the critical exponent 7 are plotted in Fig. 9. 
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Figure 9. Effective critical exponent 7 as function of the dilution p from Baker- 
Hunter analyses. 



They show an effective exponent monotonically increasing with p but reaching a 
plateau at 7 = 1 for dilutions between p = 0.42 and p = 0.46. The following sharp 
increase is to be interpreted as due to the crossover to the percolation fixed point 
at p c ~ 0.75, T c = 0, where a x ~ ex P(V^) behaviour is expected. 

It is well known (see, e.g., Ref. [40]) that series analysis in crossover situations 
is extremely difficult. If the parameter p interpolates between regions governed by 
different fixed points, the exponent obtained from a finite number of terms of a 
series expansion must cross somehow between its universal values, and does this 
usually quite slowly. Therefore it does not come as a surprise that the Monte Carlo 
simulations quoted above see the onset of a second order phase transition already for 
smaller values of the disorder strength p. The mere existence of a plateau in j e s(p), 
however, is an indication that here truly critical behaviour is seen. It is governed by 
a fixed point for which we obtain 7 = 1.00(3). Here, as always in series analyses, the 
error estimates the scattering of different approximants. 

5. Discussion 

We have implemented a comprehensive toolbox for generating and enumerating 
star graphs as required for high-temperature series expansions of quenched, disor- 
dered systems. Monte Carlo simulations of systems with quenched disorder require 
an enormous amount of computing time because many realizations have to be simu- 
lated for the quenched average. For this reason it is hardly possible to scan a whole 
parameter range. Using high-temperature series expansions, on the other hand, one 
can obtain this average exactly. Since the relevant parameters (degree of disorder p, 
spatial dimension d, number of states q, etc.) can be kept as symbolic variables, the 
number of potential applications is very large. 

Here we presented analyses of the susceptibility series for the three-dimensional 
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bond-diluted Ising and 4-state Potts model. The resulting phase diagrams in the 
p-T-plane are in very good agreement with recent Monte Carlo results. As far as the 
critical exponent 7 is concerned, however, large crossover effects render a reliable 
determination from series expansions up to order 19 respectively 18 very difficult. 
In the Ising case we estimate values that are clearly different from the pure case but 
exhibit a pronounced dependence on the degree of dilution. For the 4-state Potts 
model with its strong first-order phase transition in the pure case, the singularity 
structure of the disordered model is even more involved. Still, by comparing the 
series expansions with numerical data we can identify signals for the onset of a 
softening to a second-order transition at a finite disorder strength. 
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